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ABSTRACT 

The combined phenomena of forced convective motion 
and of transient burning of a liqtiid fuel droplet in an 
oxidizing atmosphere is studied, ffiie analysis is restricted 

f 

to the stagnation point governed by the imsteady boundary 
layer flow equations which are solved numerically to obtain 
detailed flow field for pentaborane droplets. The effect of 
various parameters on the unsteady burning process is 
determined It is found that the non-equilibrim vaporization 
kinetics plays the most dominant part in determining the 
transient behaviour. Some discussion of the time scales of 
different physical processes present in this problem is also 
given. 
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CHAPTER I 


INTROBUGTION 

I.A. SEWRRAL BAGEGROUHI) 

The study of the combustion of liquid fuels plays 
a -very important role in the design of several engineering 
systems like industrial furnaces, diesel engines, jet engines, 
liquid propellant rocket motors, etc. In fact the liquid 
fuels form a considerable portion of the total energy 
requirements of the -world. Therefore an understanding of the 
combustion phenomena becomes essential for efficient as well 
as pollution-free working of these devices , 

As an illustration of the motivation to investigate 
the combustion process of liquid fuels, the phenomena of 
combustion instability is described next very briefly. 
Liquid-propellant rocket eiagih## of various types are generally 
designed to deliver a fairly constant steady thrust for a 
duration that may extend from several seconds to as much as 
few hours. However, in many cases the resixlt is not the desired 
steady operation and instead combustion often becomes rough 
■under different circumstsmces, Ihis is mainly due to large 
pressure oscillations in the combustion chamber. Frequencies 
ranging from 10 cycles/sec to as much as 5000 cycles/sec 
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have been observed in the osoillographs. In most of the cases 
these sustained, oscillations were not truely periodic in 
nature over a wide frequency range, but seemed to be merely a 
series of random fluctuations. Not only is the thrust 
vibration, which restilts from the aperiodic oscillation in the 
combustion chamber, undesirable from the standpoint of possible 
damage to the structural elements or instruments in the vehicle 
but in extreme cases failure of the powerplant itself can occur. 
In the a^ovQ-menti^^ned situation^ the uj^gfoady 
operation is attributed to tt^e so called combustion instability, 
which is a complex phonomenoi^ involving titio inijer^OtlOP- of bhe 
combustion process with flow oscillations within the cpmbustor, 
Combustion instability initiated by a flow disturbance res'ults 
in a perturbation of steady state combustion process. Conse- 
quently the rate of energy released by the combustion process 
may increase and a part of it is fed back into the initial flow 
disturbance , resulting in a form of coupling between chamber 
pressure and energy addition to the chamber under a proper phase. 
The feed back loop so established between the combustion process 
and the flow oscillations results in a self-sustained oscillation 
of the gases within the combxxstor. In the majority of the 
cases the occurrence of combustion instability results in the 
destruction of the engine or failure of the mission supported 
by the engine. Therefore it becomes essential to investigate 
the detailed mechanism of combustion energy release and its 
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coupling to the accoustic phenomena in order to suppress this 
instability. 

In most of the propulsion devices, the liquid fuel 
is introduced into the combustion chamber in the form of 
liquid droplets of varying sizes. The treatment of the 
combustion of liquid sprays is a very complicated problem 
involving several poorly known variables like distribution of 
the sizes of the droplets, density of droplets in the spray, 
etc. Ib has been found convenient both theoretically and 
experimentally to investigate in details the burning of an 
isolated droplet, and then to apply these results for 
obtaining the burning characteristics of liqtiid sprays. 

Earliest investigations dealing with the liquid 
droplets concerned with the vaporization rates only (l). 

However during the past two decades the problem of liquid 
droplet combustion has been analyzed 'extensively under different 
sets of simplifying ass-umptions (see, for example, references 
2-14). A brief description of the physical phenomena and the 
assumptions frequently made in -determining the steady-state 
burning rate of a liquid droplet in a stagnant atmosphere will 
help to clarify the nature of the problem (8). 

A typical theoretical analysis is based on the 
following assumptions: (i) Spherical symEietry; (ii) Ho body 
forces; (iii) Pick’s law of binary diffusion, (iv) Negligible 
radiant heat transfer; and (v) The transport parameters 
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depending upon various transport coefficients as -well as the 
specific heat of the mixture are considered to be independent 
of both temperature and composition. Suitable constant values 
are often used for these parameters. 

A schematic model of the burning droplet is shown 
in Fig. 1. The physical processes resulting in the droplet 
combustion can be briefly stated as follows. The gaseous 
vapours move away from the droplet surface, mix with the 
gases of the surrounding environment via diffusion and 
convection, and then result in a zone of intense chemical 
activity (flame) at some distance from the liq.uid surface, 

A steady state is established in which a continuous supply of 
fuel (or oxidizer) vapours from the liquid stsrf ace , and the 
oxidant (or fuel) from the outer boundary sustain the 
exothermic chemical activity in the flame zone. Heat is 
transferred from the flame zone to the droplet surface to 
supply latent heat of vaporization of the liquid fuel enstiring 
a continuous supply of fuel vapours. A schematic represen- 
tation of the detailed temperature and species mass fraction 
profiles is given in Fig. 2. 

Early experimental work carried out by a number 
of workers (see, for example, references 3 and 4) had indicated 
that the relation between the radius of the droplet and time 
was well represented throTighout the steady state combustion 
of the droplet by the so called 'd law'. 


According to this 
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law, the time rate of change of the square of the droplet 
diameter remains constant (known as the evaporation constant) 
which means that the mass "burning rate is directly proportional 
to the droplet radius. TJtider the assumption of infinitely 
fast reaction rate (also known as flame sheet approximation), 
the flame zone reduces to an infinitesimal thickness and is 
represented by a surface. As a consequence of this approxi- 
mation, some explicit expressions have been obtained for the 
burning rate of fuel droplets and the temperature at the 
flame front uhder the assumption of constant as well- as 
variable transport properties (5, 6, and 14). A typical and 
widely quoted solution, based on temperature independent 
transport properties is given below (8), 

m s= "pD jtn (1 + B) ... (I-l) 

where B is the so called transfer number given by 

® = k ••• 

Ihe other symbols are defined under the heading Nomenclature. 
These theoretical results agree very well with the experimental 
observations bn steady-state burning.. 

However, almost in all practical spray combustion 
systems, the biuning takes place under the forced convective 
conditions, and in addition, the combustion phenomena, for an 
appreciable fraction of the total droplet lifetime, remains 
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transient. It has been found that the pollutant formation 
in auto-engines, and the comhustion instability phenomena are 
generally affected significantly by the transient processes. 

In the present analysis an unsteady model was used to treat 
the biorning of liquid fuel droplets in an oxidizing atmosphere 
under forced convection. Before going into the details of 
the present analysis, it may be xiseful to review some of the 
important past contributions dealing with these two aspects 
of the droplet combustion, 

I.B, REVIEW OP IHE RELATED PAST aOITTRIBUTIOIiS 

Although a number of theoretical analyses have been 
done on the combustion of isolated single liquid fuel droplets 
in an oxidizir^ atmosphere during the past two decades, there 
are still very few attempts dealing with the combined phenomena 
of transient burning and of forced convective motion. A brief 
description of the past work investigating these effects 
separately is given below. 

At the very early stages, it was realized that 
the spherical symmetry of the droplets could not be achieved 
because of the presence of natural or forced convection effects. 
An early attempt to determine the influence of forced convection 
experimentally was made by Grohrbandt (9). The mass evapora- 
tion rate of camphor spheres was found proportional to the 
radius and the square root of Reynolds n-umber. Spalding (3) 
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applied the concept of "boundary layer to take into account 
the Gontri"butions of both natural and forced convection. He 
established an empirical relation to predict mass burning 
rate under the influence of forced convection by performing 
experiments on spheres over a range of Reynolds number varying 
from 800 to 4000. Several other empirical relations incor- 
porating forced convection have been found out since then (14). 

Ror example, more recently Paeth and Lazar (12) obtained the 
following correlation in terms of the burning coefficients, 

^(forced donveotlon) * 0.278Rel/2 Pr^/^d + , 

... ( 1 - 5 ) 

tdiich was found consistent with the experimental data in the 
range 10 Re < 800, and also the theoretical prediction of 
Pendall et.al (lO). Recently Bruzustowski and Hatarajan (ll) 
studied experimentally the evaporation and combustion of cold 
n7=spentane d^’eplbts injected into a hot oxidizing gas at 
moderate presstires with hi^ speed motion picture photography 
and were able to correlate their results in terms of transfer 
number, Reynolds number and Husselt number. 

Simultaneously there have been some attempts to 
study separately the transient aspects of the problem. 

Several of them deal with Idle supercritical conditions of the 
liquiid phase such that the droplet is treated simply as a 
puff of gaseous material diffusing into the surrouinding 
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enTiroiment (14-19). Dae 'b-urning of droplets under such 
conditions is highly transient in nature and moreover the visual 
theories to determine the combustion phenomena becomes highly 
inaccurate. Spalding (15) first considered this case by 
treating the liquid droplet as a point source of gas injected 
instantaneously into an infinite medium and predicted the 
bvirning rates Tdiich are valid at hi^ pressures, Rosner (16) 
modified Spalding’s results by considering finite size gas 
sources, later Ghervinsky (17) extended the analysis by 
including the convective terms and the variable density. He 
obtained solutions of the flow field arpvind a burning droplet 
in the Von-Mises plane using flame sheet approximation. 

Recently Polymer opoulos and Perkin (18) extended this model 
further by including the finite rate chemistry. Rosner and 
Chang (19) investigated the transient evaporation and 
combustion of a fuel droplet near its critical temperature 
and showed the inadequacies in the past analyses based on the 
familiar quasi-steady approximation. 

Ihere are a few ^vestigations dealing with other 
transient phenomena associated with liquid phase evaporation 
or combustion in stagnant enviroment under normal (non-critical) 
conditions (14,15,20-23). Rirkaldy (20) first tried to solve 
the time dependent heat and mass transfer equations by 
assuming slow evaporation so that the convective term in the 
equations disappear. His analysis was restricted to the 
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ph-enomena of condensation without any chemical activity. 

Spalding (15) also obtained analytical expressions for the 
combustion of a fuel droplet at high pressures by assuming 
xmiform properties including constant density. Strahle (21) 
employed perturbation method to the unsteady spherical symmetric 
burning problem by representing the dependent variables as 
power series in -vdiich the leading terms were time independent. 
However, the first order time dependent solution was obtained 
assuming a constant density throughout the flow field. 

Chervinsky (22) obtained the detailed solutions for the fuel 
mass concentration and temperature fields in the Von-Mises 
plane as well as the burning rate as a fimction of time by 
employing flame sheet approximation. Kotake and Okazaki (25) 
obtained nimierical solutions for ^lnsteady vaporization and 
combustion of a fuel droplet in quiescent air assuming infinite 
kinetic rate reaction. It is shown that the process of 
evaporation and combustion are influenced greatly by the 
initial conditions of the air surroxuading the droplet and 
that the combustion process approaches a quasi-steady state 
near its end while the process of evaporation cannot attain 
a quasi-steady state of mass and heat transfer at the droplet 
surface. 

To investigate the combined phenomena of transient 
burning and of forced convective motion of the oxidizing gas 
relative to the droplet, Strahle (24) employed the unsteady 
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■boundary layer flow approximation to the forward stagnation 
point of a Vaporizing droplet and then obtained approximate 
solutions by using an integral approach, like Eotake and 
Okazaki, it was found that the development of thermal layer 
within the liquid did not attain a steady state dxiring the 
typical droplet lifetimes. For aecurate analyses, the 
transient vaporization kinetics must "be included. 

In the present analysis a ftilly unsteady model is 
used to deal with the oombined phenomena of forced convective 
motion and of transient burnii^ of a liquid droplet. The 
treatment fs however restricted to the stagnation point in 
view of the resulting computational simplification as well as 
the fact that it can also represent qualitatively the essential 
features of the entire droplet combustion process. However, 
the procedure is easily extendable to other regions of the 
flow field aroTuid the droplet. The unsteady conservation 
equations of species concentration, velooity, and temperature 
of gaseous mixture surround|.ng the droplet and of temperature 
within the droplet ^e solved by finite difference scheme which 
is more exact compared to that of an integral approach. 

Accurate detailed profiles for velocity, temperature, and 
species concentration are obtained as a result of this analysis. 
The various time scales which play an important role during the 
transient combustion of a liquid fuel droplet are described. 

The details of the present analysis is given in the next chapter. 
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CHAPTER II 

EORMIXLATIOH OP THE PROBLEM 

II. A. PHYSICAL MODEL 

In general, the liquid fuels or propellants are 
Injected through orifices into the combustion chamber in the 
form of liquid jets. These jets soon breakup into tiny 
droplets. Upon ignition of the spray, droplets vaporize 
resulting in a supply of gaseous fuel and oxidizer. These 
vapours mix with each other and react forming hot combustion 
products. Consider a single fuel droplet in the spray which 
undergoes transient burning in a forced convective oxidizing 
environment of the combustion chamber. Hue to the flow of 
hot oxidizing gases at temperature T. and pressure P with a 
velocity U^ relative to the droplet, a boundary layer is 
developed OTner jjhc auffaee of th® liquid fuel droplet (see 
5)> However, in the wake region, the boimdary layer 
approximation breaks down and the analysis of the flow field 
becomes highly complicated. The boundary layer flow portion 
is treated as laminar for the purpose of this analysis. In 
reality the turbulent conditions are more likely to prevail. 
However, this can be taken into account phenomenologically 
in the numerical analysis developed for the present model by 
the introduction of proper exchange coefficients for viscosity, 
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thermal conductiTity , and coefficient of diffusion in the 
governing equations , 

Initially, after ignition, the fuel evaporates from 
the liquid stirfase, diffuses outwards and mixes with the 
gaseous oxidizer diffusing mwards from the surrounding free 
stream. With the passage of time heat is liberated as a 
result of exothermic reaction taking place between the gaseous 
oxidizer and the fuel vapour inside the boundary layer. The 
heat so liberated further augments the process of vaporization. 
Therefore the transient burning of a liquid fuel droplet in a 
forced convective environment involves the solution of tmsteady 
boundary layer flow equations coupled with the heat and mass 
transfer conditions at the surface, I’or simplicity, the 
present analysis is restricted to stagnation point of the 
liquid fuel droplet. 

As far as the droplet interior is concerned, it is 
assiimed that no circiilation of fluid occurs and that there is 
no mixing of the liquid fuel with the oxidizer or product 
species. In other words, only heat conducted within the fuel 
droplet is included in the aioalysis. Moreover, it is also 
assumed that at the centre of the droplet the temperature 
does not change >±Lich implies that the thermal diffusion layer 
at the liquid fuel droplet surface is much thinner compared 
to its radius. Also, when the analysis is restricted to 
the stagnation point, the temperatiare gradient along the 



15 


droplet surface is neglected compared to the gradient in the 
perpendicular direction, 

Ihe flow configuration representing schematically 
unsteady burning process analysed in the present model is 
shown in Fig. 4. The origin is fixed at the stagnation point 
0. X and y axes represent the distances along and perpendicular 
to the droplet STorface respectively. Dotted line shown in 
Fig. 4 indicates the edge of the boundary layer. It is clear 
that the stagnation point flow approximation reduces the 
spatial dependence of the partial differential equations by 
unity^ Due to ssmimetry, all the dependent variables can be 
chosen to be functions of y alone, except U which must be 
taken proportional to x. 

The detailed mathematical formulation is described 
in the next section. 

II. B. MATHEMA.TIG1D FORMULATION 

II.B.l SQVprm' ng Equations ? 

The governing equations for the chemically reacting 
two— dimensional, unsteady boundary-layer flow can be expressed 
as follows (see, for example , ref erenee 8) ; 

Continuity ; 

+ pv.u = 0 II-l 



16 





i=l "i 

And the total deri-vative. 


Ei_i _ . riLl+ „ lJLi+ ^ JLl] 

Dt “ Ljt » sy •' 

The symhols used in the above equations are defined under the 
heading Nomenclature. 

The continuity equation is simplified by neglecting 
the local variation of the mixture density with respect to- 
time as compared to its spatial gradients. Then Bq. (II-l) 
reduces to 

iL Cp^) 4. L- (p v) ^ Q 
«x 3y 

Since the pressure gradient in a boundary layer 
is impressed by the external stream, we have 
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dU . 

p jj ^ ^ 

^ "W- ^ >w 


{ 11 - 8 } 


e e^dx ■' ~ ~ dx 

Substituting Bq., (II-8) in Eq, (II-2), the momentxmi conservation 
equation becomes 


,£il 

Dt 


dU a ^ 

_ p u ( ) + r~“ ( u, 
e e dx ^ 9y 


(11-9 ) 


In general expressions for the y-comp6nents of the 
heat flux vector (q), and of the diffusion velocity of species 
i (7^) are knoun from the molecular theory of gases in terms 
of the other dependent variables. Here we will consider highly 
simplified form of these expressions as discussed below. 

Oonsider the general expression for diffusion 
velocity for a multicomponent gaseous mixture due to concen- 
tration gradients only (8) 


H 

! 

j=i 


¥. D. . 


(V, - 7,} 
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H Y- VY 
J -1 f—J. 
4 w ^ Y 
0=1 0 ^0 


V Y. 

i 

^i 


) 


i 7^ o' > i = 1,2,. . . 


(II-IO) 


If we assme all and to be equal, then the above 
equation reduces to 


Y Y 
^i''i 


- D. YY. 

X 1 


(II-ll) 


which is the well known Pick’s law for binary diffusion. 
Therefore the y-component of diffusion velocity in Bq. (II-4) 
is expressed as 


D 


8Y. 


7. = - f- (t^) 


(11-12) 
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Hext, the heat fltix vector for a multicomponent 
mixtiire comprises of two mechanisms, namely, the heat 
conduction and the thermal energy carried by inter diffusing 
species. This can be expressed as 

^ ^ 

q = - X V T + p I h. Y. t. (11-15) 

i=l 1 m 1 

Here we have neglected the contribution due to thermal 
diffusion of species. Therefore, the y-component of the heat 
flux vector, q, in Eq. (II-3) is approximated by 

^ ''ill 

where is given in Eq. (Il-rlS). 

It is convenient to express the energy conservation 

equation, Eq. (II-3), iu tems of temperature (T), and for 

this purpose , it is further assumed that the specific heats 

at constant pressure of different species are nearly the same 

and can be replaced by an average value independent of 

temperature, that is, G = C = constant? for all i. Eien 

P P 

making use of Eqs, (11-12) and (11-14), we can rewrite Eq. 

(II-3) as follows 
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In order to express the soinrce terms involving 
in Bqs. (11-4) and ( 11-15), we approximate the chemical 
reactions taking place in the flame zone by an overall single 
step reaction of the following type: 


V®] + '’I'M * + ••• (11-16) 

fhe order of the reaction for a given combination 

of fuel and oxidizer is determined from the kinetic information 

quoted in the literature. Now, according to the law of mass 

"felx 

action, the net ra,te of production per unit volume of the i^ 
chemical species | w^, is given by 


N 


V . 


w. 


n t . Ml" ''A 

= W.(v. - V )w = ¥,(v - V, )k n (c.) J 


D=1 


i = [0], [F], (11-17) 

where the specific reaction rate coefficient k is approximated 
by the Arrhenius expression, 

k = (11-18) 


H 


we will let a suitable effective value ^ 


3=1 


+ P-p) 


to represent the order of the reaction in the following analysis. 

Having obtained a detailed expression for w- , we 

N : : ^ V 

can now express the term | ^i'^i terms of other variables 


i=l 


zs follows: 


N N , 

I = wj W,(v, 


i=l 


i=l 


)h^ = w q^ 


(11-19) 
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Here q.^ stands for the heat of comhustion at the reference 
conditions. In Tiew of Eq.s. (11-12), (11-17), (11-18), and 
(11-19), the species conservation equations, Bq, (II-4) and 
the energy conservation equation, Eq. (11-15), become 


Ho 


3 


- ( 


w. 


_) 0""° /° „-SAl°T 


'’0 ^ 


Tj, e 


'0 


w. 


(11-20) 


3Y. 
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( 11-21 ) 
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^ 


( 11 - 22 ) 


It is convenient and consistent with the simplifi- 
cations introduced above to assume a constant average molecular 
weight (¥) for the mixture, Then, the equation of state for 
ideal mixture, Bq, (II-5), becomes 

p = p fi®T/¥ (11-23) 

Next equation governing the heat transfer within 
the liquid fuel droplet is discussed. Neglecting the tempera- 
ture gradient in the x-direction compared to that in the y- 
direction, the unsteady heat conduction equation for the 



liquid droplet can be expressed as 


3t “ ^ ^ "ay”^ (11-24) 

Fote that the droplet surface is regressing inwards at a 
rate equal to the vaporization rate. Per the non-equilibriuia 
vaporization rate, the mass lost by evaporation per unit time 
per unit surface area can be expressed as (8) 


Ps^s = CP 


""h s^ ' s3 

27tR°Tg eq ^ 


(11-25) 


here (l-rp _) denotes the equilibrium mass fraction of the 
^ eq 

fuel vapour at the surface corresponding to the temperature 
and pressure p. £ , the evaporation coefficient is defined as 
the ratio of the observed rate of evaporation to the evaporation 
rate calculated from the Khudsen equation. When the liquid is 
in equilibrium with the vapour, the equilibrium vapour pressure 
of the fuel (Pp g^) ^ related to the temperature of the 

liquid (2^) with the help of QlSpeyron equation as follows (8) 


■^o ' 

^5',eq 


= exp [ 


(1 - 


( 11 - 26 ) 


where pS _ denotes the value of p-c, _ at the initial temperature 
(!I^ q) of the liquid. Fote that the above relation is applicable 
when both phases contain only one species. Purtheimore , it is 
assumed that is constant. In otherwords, any temperature 
dependence of the heat of transition per mole (H ) has been 
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neglected. Thus (Yp in Bq.. (11-25) can be expressed 

* eq 

with the help of Bq. (11-26) as follows: 




eq 


(Yp ) exp 
^ eq 


Eo 


R®T 


(1 — 


T 


fctO 


)] 


( 11-27) 


5..0 


0 

where (Yp denotes the equilibria mass fraction of fuel 

^ eq 

vapours at the surface corresponding to the initial temperature 
of the droplet T^. 

In order to fix gas-liquid interface at y = 0 the 
liquid surface is continuously translated towards the interface 
at a velocity equal to the instantaneous regression rate due 
to vaporization (v ). Due to incompressible nature of the 

£ f S 

liquid any interface pressure fluctiiations caused by this 
motion are considered negligible . Introducing a modified 
;’'-coordinate , the unsteady heat conduction equation within 
the droplet, Bq, (11-24), then becomes 


3T£ STj, 3 aT, 

[— ■ 

t > 0, 0 < y^ < - » (11-28) 

,t 

where y « y + / v„ „ dt. [Note that y (or y ) is negative 
S. Q SLfS * 

towards the liquid phase side.} At the surface, from mass 
conservation, we have the relationship 


(11-29) 



II. B. 2 Boundary Conditions; 


The boundary conditions for continuity and momentum 
equations are easily expressible in terms of the velocity 
components as follows; 

At the droplet surface, y = 0; 

u s! 0 ( ITo slip condition) ( 11-30 ) 

V = v^(x,t) (Due to vaporization) (11-31) 

o 

And, for away from the liquid droplet surface, y-». « : 

u = u(x,t) (II-32) 

The boundary conditions for the chemical species 

at the gas liquid interface (y = 0) are determined by mass 

balance at the interface. In general, the component i is being 

transported towards the gas-liquid interface by diffusion at 
8T. 

the rate (pD^ ’ and, simultaneously it is being converted 

s 

away from the interface by a normal current at the rate 
(pv) (Z- ) and towardv the iuterface at the rate (pv) (Y. ) 

provided, in the last case, the component i exists in the 
liquid phase. Thus, the boundary conditions for the oxygen 
and fuel mass fractions are respectively given by 

(11-33) 

S' ■ 
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where we have ass-umed in Eq.. ( 11-34) that the liquid droplet 
contains only pirre fuel. Outside the boundary layer (y — »- » ), 
it is clear that the species mass fractions must satisfy the 
following conditions: 

Yj, = 0 , and Yq , Y^^^ (11-35) 

The boundary condition at the gas-liquid interface 
for temperature is determined by considering energy balance at 
the Interface. Ihe heat conducted from gases to liquid surface 
is equal to the sum of the heat conducted from liquid surface 
to the interior of the liquid and the heat absorbed for 
vaporization at the surface 

9T 

l(-x^) I = I (- X, I + H, (11-36) 

o o 

Again, as y — ^ ~ , we have 

T = T (11-37) 

The liquid fuel droplet near the stagnation point is assumed 
to be like a semi -inf inite mass without any internal circulation 
(see Fig, 4) . Therefore, the temperature of the liquid far 
away from the surface (yj^ remain at its initial constant 

value, that is, . 

( 11 - 38 ) 

This approximation is quite reasonable as long as the thermal 
wave penetration thickness is much smaller than the droplet 


radius , 
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The fuel droplet is instantaneously introduced into 
the oxidizing environment. This represents a discontinuity 
in the physical variables like velocity, species mass fractions 
and temperature which req.\iires special care in specifying the 
initial conditions in order to start the numerical integration, 
and is described later on. 

II. 0 TRAHSPORI^IITION OF GOYERNIITG EQUATIONS 

It is convenient to solve the boundary layer 
equations in terms of the well known independent variables 
introduced by Howrath and Dorodnitzyr as well as Levy and 
lees in the following form (See, for example, reference, 25) 

X 

5(x) = / dx (11-39) 

0 

u J 

Ti(x,y) = -n=. / P dy (11-40) 

m 0 

In addition by introducing a stream function defined by 
the relations 

PU = If . a.a PV = (11-41) 

the continuity equation, Eq. (11-7), is identically satisfied. 
However, it is common to use another modified non-dimensional 
form of the stream function, namely, 

f(c,n,t) = ^ (11-42) 

m-' 
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Then, from Bq.s, (11-39) to (11-41), it can "be easily shown 
that 

(^) = H (11-43) 

© 

XT 

and pv « (11-44) 

The time variable (t) is also non-dime ns ionalized by dividing 
it with a characteristia flow time (t^)> 

i: = t/t^ (11-45) 

The choice for t^ is given below in Eq. (11-54). ITow trans- 
forming Eq. (II-9) and Eqs, (II-20) to (11-22) with the help 
of Eqs. (11-39), (11-40), and Eqs. (11-42) to (11-45), the 
conservation equation for momentum, fuel and oxidizer mass 
fractions, and temperature, respectively, become 




o^ ri£i-r^\ M iff i 

^ ‘-STl 35 '‘3^ 31 3^2-' 


( 11-46 ) 


+ -P „ i-i 

S-n^Sc ^r]‘ BT) ^ 3v 


Bp 4 25 [Cf - f -^1 


3 jHo.) ^ f jl!^ !5 q. _ b + 2e r^..^ ^ 


] ( 11 - 48 ) 


3 fC ad'v , « , e/ p jL Y 

sn^Pr 3T)' ^ 3n 911^^ 




3f 3© 3f 3©, 


- - QBp + 25 ^ 3X1 *35 "" 35*3 ti 


(11-49) 
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Here C = ( ■ ) , Pr >= , and Sc = ([i/pD) are, respectiTely , 

^e^e 

known as the Chapman-Ruhesin function, the Prandtl number, and 
the Schmidt n-umher, Q =; Cq. 3 ,/( 3 represents the non- 

dimensional heat of combustion. The source (or sink) terms 
Bp and Bq are given by 




A t, 
a 




f-,S- 1 /o /p ^-e/© 

P -1 J ^ ^0 ® 


*0 *P 


n 1-1 


0 


A P 1*— n (X p-CT ^ /q 

V r ^ S— 1 a 0 Y ^ Y ^ 

^0'- a_-l Pp •* 0 ® 


V ^P 


( 11 - 50 ) 


(11-51) 


Por the heat transfer equation within the fuel droplet, Bq, 
(11-28), we introduce another dimensionless variable in the 
y-direction. 


“n, 


= I 


PoG 






( 11 - 52 ) 


So that Eq,(II-2Q) becomes 


3T) 


2 


- V. _ V — ; ' 


l,S 




= 0 


( 11-53) 


t * 


where we have assumed a constant > 

In the region of the stagnation point, the governing 
equations can be highly simplified by dropping 5-dependence 
altogether. Moreover, the velocity field for the external 
flow can then be expressed as 


Ug = (x/t^) , and t^ “ st.pt = constant (11-54) 
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Now the governing Eqs. ( 11-46) to (11-49) reduce, respectively, 
to the following forms: 


£_ y.' 

Sc 


y,. 

■Sc ^0- 
0 


fp' + 

(© - p^) 

- = 0 

St 

(11-55) 

+ fYJ, 


=1- 

(11-56) 

+ fY- 

.ik _ 

9T 


(11-57) 

■ f©’ - 

3t "" 

QBj, 

(11-58) 


Here ( )' denotes differentiation with respect to r). Furthermore, 
the expression for the mass flux of the fuel vapours from the 
droplet surface Eq.. (11-44), reduces to 


^s s 


^e ^e ^ __ / Pe ^e * p 

- “ f t^ s - - / C_ t^ ^s 


(11-59) 


.f w , -.p vy 

where, for the last equality sign, Prandtl number is taken as 
unity. Now, making use of Bqs. (11-31) and (11-59), Sq. ( 11-53 )■ 
becomes 


— = ° • ^>0, 0< y, < 


9 N 


Cp 9= 


(11-60) 


where 






( 11-61 ) 


Note that y is negative in the above equation. Similarly 
making use of Eqs. (11-31) and ( 11-59), Bq* (11-25) can be 


expressed as 
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) ( 11 - 62 ) 

Rewriting Eq. (11-27) in terms of non-dimensional symlDols, 
we obtain 

= (Yj, )° eKp[e^(l - (11-63) 

' eq * eq 

Having transformed the governing eq-uations, the 



boundary conditions should also be expressed in terms of the 
new variables. Following the same procedure, Eq. (11-30) to 
Bq. (11-37) then become. 


At the droplet surface, r] = Oi 



(11-64) 

(11-65) 

(11-66) 

(11-67) 

( 11 - 68 ) 

# ' ■■ 

(11-69) 

(11-70) 

(11-71) 


1 


(11-72) 
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Also, the houndary condition in the liq.'aid phase, Eq. (11-38), 
can he rewritten as 

At ,1^ — : S, = (11-73) 

Eor the initial conditions, the solutions obtained 
in reference (24) for an impulsively produced boundary layer 
flow have been employed. This final form of the transformed 
set of governing equations are solved numerically, and the 
method is described in the next chapter. 
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CHiPTSR III 
METHOD OP SOLDTIOH 

III. A. FINITE DIPPERENGE EQUATIONS 

In order to solve Eqs. (11-55) to (11-58) and Eq, 
(11-60), which form a set of coupled nonlinear partial 
differential equations, each equation is first transformed into 
the finite difference form. Implicit method was employed 
to obtain the finite difference form of three equations. Thus, 
a spatial partial derivative is approximated by finite 
difference terms evaluated at the advanced time The 

transformation is achieved by establishing a network of grids 
throughout the region 0 t) and 0 < v £ '^xaax 

spacings an and at respectively (see Pig. 5). Pig. 6 shows 
schematically the grid points involved in the space and time 
difference approximations xased in the analysis. At any grid 
point denoted by (i,n)’, the finite difference form of the 
governing equations can be represented by the general' equation 

“■^i ^(i-l,n+l) ®i^('i,n+l)'‘ V(i+l,n+l) ^■(i,n)» 

2 < i £ M (III-l) 

At any one time level , Eq* (III-l) will be written once 
for each grid point 2 i M. This results in a system of 
(M— 1) simultaneous equations with (M-1) unknown variables in 




^(i,n+l)* Before going into the method of solution of such a 
system we will obtain the finite difference form of the 
governing equations. 

Since Bqs. (11-55) and (11-58) are nonlinear 
equations , the nonlinear terms are linearized in the following 
manner . 


Consider a nonlinear function G-(gj>g 2 » • • • »%) of 
several dependent variables like g^jg^ , . . . ,gjj. This function 
can be expanded in a Taylor's series about a value of g^ 
either guessed or knowi from previous iteration. That is. 




>%) ] 


(k) 


0=1 


’) ( 


3G (k+1) 


- 4 “^^) 


(III-2) 


where k represents iteration number for evaluation of G-. Note 
that as soon as becomes equal to g^^^ within a prescribed 

limit, and als# become equal to that accuracy. 

With the help of Eq. ( III-2) , the nonlinear terms 
of Eqs, (11-55) and (11-58) can be expressed as follows: 


.exp{-e/G^, „^-,^>_(k+l) (ij.) ^ -.(k) Q(k+l) 

(i,n+l) (i,h+l) (i,n+l) 

(II 1-4) 


.(k) 




^(i,n+l) 


= a 
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where =■ (™-5) 

ra's-'^^ 1 0 


and 


■u(k) 

^(i,n+l) 


_ i.n+l)\ r v^(k) 1 

>-‘^(i,n+l)J (111-6) 

The various derivatives in the governing equations 
have been approximated by finite difference expressions as 
described below. 

The first order derivatives with respect to h like 
p', Y^^ ©' and are approximated by the central difference 

formula. That is, 

Similarly , fur p” , YJ , Y^ , ©" and ©J , the following three 
point difference approximation is used, 

,(k+l) n. 


P_(k+i) . xl/h^ 

^(i n+l) ^^(i+l,iH-l) “ 2P(i,n+l) P(i-l,n+l)-^/^ V 


etc. 

To express the time derivatives, (ap/S'f), etc., in difference 
form, the backward difference scheme is used as given below, 

(k+l) 


(«) 


aT' 


(i,n+l) 


_ ij i P-tt ) , v(t|n)„ j ^ etc. 


(III-9) 


AT 
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i:ot9 that the q.uantities already known at the previous time- 
step do not require any iteration index. In order to evaluate 
the first order derivatives of the dependent variables at the 


droplet surface which appear in the interface conditions, 
sli^tly more accurate three point difference formulae have 
been used. That is, 

— i— f/lY _ Y „ ■:?Y (^+1) 1 

° (l,n+l) °(2,n+l) °(3,n+l) ^(l,n+l) 


(III-IO) 


( k+1 ) 

^ (l,n+l) 


^ 2h l-^^(2,n+l) ^(3,n+l) ^^{l,n+l)^ 


2h 


[4Y. 


(k+1) 


- Y 


(k+1) 


— — 3 Y ^ ^ 1 

^(2,n+l) "^(3,n+l) ^(l,n+l) 


(Ill-ll) 


(l,n+l) 


(II 1-12) 


(k+1) 


.(k+1) 


^ *ll,n+l) ^ ^ t2.n+l) “( 3 , 11 + 1 ) ' "(l,n+l) 


_ _ ^gVx+.-^Z J 


(III-13) 

Therefore Iqs. (11-55) to (11-58), and Eq. (11-60) are now 
readily transformed into the following set of linear finite 
difference equations: 

,2 -(k) iJt+l), 


- r[C - ^ [1 + src + 2rh2 

- + (V2)f^f,„+l)]p[lrt!n+l) 
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-rCC/Sc - (h/2) 


(k+1) 

(i-l,n+l) 


, (k) 




+ [1 + 2rC/Sc + 


& 


( i ,n+l) 


1 y (k+l) 

I ^F(i,n+1) 


-r[0/Se + (h/2) Tj, - 


(k+l) ^ Y 

(i+l,n+l) ^(i,n) 


(111-15) 


-r[0/So - h/2 ^0 


(k+l) 


+ [l + 2rG/Sc + 


( i-l,n+l) 

Vo¥orh"K(Tj , 


2-'- ) exp{-«/4?^ 


^k) 


r\\^. 

*'(i,n+l) 


, „ (k+l) 

°C 1,0+1} 


-r[0/So + (h/2) To 


(k+l) 

(i+l,n+l) 


= I, 


0 


(i,n) 


( III-16) 


-r[0/Pr - h/2 


,(k) 


+ [1 + 2rC/Pr - VjWjrh-KQ 


(k+l) 


-r[0/Pr + (h/2) 


— ^ rh-^KQ ^ >^,4.1 'i ^F , 

(i,n+l) 0(i^n+l) ^(i,n+l) 


( III-17) 



37 


c 

-r[l + f (h/2) — + [1 + 2r] 6, 

® “(1-1, n+1) »(i, 


n+1) 


-r[l - fgdi/Z) ©, 

'P 


c 9 ^ '“(i+l,n+l) 


= ©n 


(i,n) 


(III-18) 


where 


r = ^ and K ^ (At^/ p-|Wq¥j,) 


( III-19) 


¥e have assmed here an over-all single step second order 
chemical reaction between the fuel and the oxidizer vapours. 

The interface conditions stated in Sqs. (11-66) to 
(11-68) with the help of Bqs, (III-IO) to (III-13) become 
At the droplet surface T) = 0: 


P(l,n+1) 


= 0 


Lp 


( k+1) 

( 1 ,n+l) 


= Y. 
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(3 - 2hf J 
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(III-2G) 


(I 11-21) 


(III-22) 


n(k+l) _ 

^(l,n+l) “ % 


4rQ(k+l) „iQ^(k+l) ^ ]+[2hQ ] 
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(III-23) 
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And, at the edge of the boimdaxy layer, t ] — « ; 

^(M+l,n+l) ” ^ (III-24) 

The numerical procedure of solving the finite 
difference equations, Eqs» (III-14) to (III-24) is described 
in the next section. 

1114 SQXUTIQM TO THE FIHITB DIPFEREhOE EQUATIOITS 

As already mentioned each of the finite difference 
equations, namely Eqs. (111-14) to (l|lr-i8) can be represented 
in the general form of Bq. (III~1). As a result (M-l) linear 
equations will be obtained at each time step. These equations 
can be represented in the form of square tridiagonal matrix of 
coefficients A^, and C^. The right hand side of these 
equations, lj_is a known quantity. Therefore, with the help of 
appropriate boundary conditions, solutions to these tridiagonal 
system of equations can be readily obtained by (Jaussian 
elimination technique ( 26 ) . Th^s tqchi^quie is computationally 
stable for a vast choice of ^tep si^es and has outweighed the 
moderate increase in computational effort when compared to 
the qther techniques. Numerical computations were done by 
programming in Bortran on IBM 7044/1401 digital computers. The 
following important points must be kept in mind before 
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proceeding for computation. 

As already mentioned in Section III. A, the implicit 
difference scheme involves the computation of values at an 
advanced time step, making use of the evaluated values 

at the previous time step Since the droplet is injected 

instantaneously, there will he sharp gradients in the flow 
variahles y = Q, Hcwever, the ntperical treatment is made 
possihlf hy starting integration after § veny small interval 
of time I ihe initial values for q.uantities like p, Yp, T^, 

9’^ and 9 are calCi'lated from the fimcticns developed for the 
Integral technique in reference 24, A reasonable guess for x 
(t = 10“^) was made to startwith%- Initial guess was made for 
mass fraction of species (Yq ^ and Ip g)» and f^ at the 
interface*. 

Stable solutions were found out for each of the 
profiles f , p, Yq , Yp and ft, by solving the finite differenc 

equations simitltaneously . In obtaining the solution the wall 
values are assumed at first hnd then checked with the botuadary 
conditions expressed by three point difference formulae. For 
example, Bq, (111-^23) was used to obtain refined interface 

temperature (ft j* llhe process is. repeated till all the 

S', ' ' , ■ [ 

interface values converge to the desired accuracy. As soon as 
the new profiles were obtained old ones were replaced by them 
in order to compute for the next time step. The iteration 
procedure was repeated until a prescribed convergence of 
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1 in 1000 was obtained f'or all tlie netjuired qiiantities. 

In order to advance the solution in time , a sman 
increment in the time step ( at) was given. The solutions 
obtained at the previous time step were used as the input 
profiles for the next time step. The same procedure is 
repeated as described in the previous section until a 
converged solution is obtained at this T-values. The 
computations are stopped when the field variables reach steady 
state. It was found convenient to use a variable time. step. 
Time steps of 10 to 10 were used for the unsteady part of 
the solution. 

As3rmptQtic approach of the profiles p, Yq, Y^, © 
and was obtained by trying various values ef . An 
optimum choice of ® found adequate to satisfy 

the external flow conditions. 
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CHAPTER lY 

RESULTS, DISGUSSIOES AND CONCLUSIONS 

IV, A RESULTS AND DISCUSSIONS 

The computations for transient comhustion of a fuel 
under forced convective conditions was undertaken to determine 
the time scales characterizing the various physical processes 
involved in it. The numerical values of the physical constants : 
appearing in the governing equations are chosen to correspond to; 
those of stable liquid pentaborane (B^Hg), As compared to 
many hydrocarbon fuels, pentaborane is chemically more active, 
has a high calorific value and is stable over a larger period 
during storage { 27 )* Por poorly known values of some physical 
constants, a set of different values were tried. In addition, 
the values of certain parameters were varied in order to study 
their effect on the combustion process. A list of the data 
used in the numerical calcuilations is given in Appendix I. 

Bach case there refers to a set of input data for one complete 
solution. 

Detailed profiles of velocity (p’)> species 
mass-fractions (Iq and T^), and temperat\ires and as a 
function of t and rj are shown respectively in Pigs. 7 to 11 
for one set of physical data. (case I of Appendix I). It is 
evident from these figures that the gaseous boundary layer 










DIFFERENT TIME INTERVALS 


FIG,7, VELOCITY PROFILES A 








FIG 9 MASS FRACTION, PROFILES OF OXIDIZER AT 
":LTlig»WfBRVALS,- : I 
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grows with passage of time tintil the steady state profiles 
are obtained. In fact the viscous boundary layer thickness 
does not change as soon as x becomes of the order of unity 
which is reasonable since the characteristic flow time is 
taken equal to t^. The velocity and temperature profiles 
(Tigs, 7 and 10) attain steady state behaviour at about 
X = 10 while the species mass fractions profiles (Pigs. 8 and 
9) do so around x = 100. In otherwords, as expected the detailed 
behaviour of the gaseous boundary layer is controlled by the 
species diffusion rates , before the gaseous phase attains a 
quasi-^steady state behaviour. As x increases, the thermal layer 
thickness in the liquid fuel goes on increasing in an asymptotic 
manner which requires more computing time with a variable 
However, using the steady-state constant values for f and 
© , it is possible to obtain the steady-state temperature 
distribution in the liquid fuel with the help of Bqs. (11-60) 
to (11.65). The result is shown by a dotted curve in Pig. 11, 

It is clear that the temperature distribution obtained from the 
unsteady equations does not attain steady-state results even 
when T = 100, unlike the processes in the gaseous phase. Thus, 
the wet bulb temperature is not achieved at the liquid fuel 
surface. It is therefore important to note that the vaporization- 
kinetics controls the unsteady combustion of the fuel. 

It is also instructive to study the behaviour of the 
variables at the gas-liqifid interface. Pig. 12 shows the 
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-variation of the vaporization rate (- f ), the interface fuel 

mass fraction and the interface temperature (Q^) with 

time (v). The interface temperature remains practically 

constant. However, it must he pointed out that even a slight 

variation in its magnitude results in a significant change in 

the vaporization rate because of the exponential dependence in 

the kinetic terms. It has been assumed in the present analyses 

that there are no fuel vapours initially in the surrounding 

environment when the droplet is injected into it. Consequently, 

the vaporization rate must be high at the initial stages such 

that the liquid provides a part of its own latent heat of 

vaporization which results in a sli^t decrease in the magnitude 

of the interface temperat-ure. The starting values of both the 

■fuel vapour pressure (Y-p „) and the evaporation rate (- f ) 

CQ, • 

are determined by the initial value of the interface temperature 
(&„) which is chosen to be between unit and ©. . Note that 

Yp g will be close to zero at the initial stage because there 
will be hardly any fuel loading into the gaseous boundary 
layer. However, Y-p starts building up as time increases 
(see Pig, 12), As r approaches mity, the boundary layer growth 
becomes steady as described previously but , from Pig, 12, we 
find that the interface temperature begins to increase -vdiich 
results in a sli^t increase in -fg and equal to 

about 50, interface temperature ceases to increase, and 
remains constant afterwards. Instead, keeps on increasing 
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till X is around 100. Thus the interface fuel mass fraction 
does not respond instantaneously to the increased vaporization 
rate, and there is certain time lag within which ^ builds 
up corresponding to the new vaporization rate. 

Results of the computations based on an integral 
method (24) are shown by dotted curves in Pig. 12. There is 
an appreciable difference in the absolute values of the different 
variables. Also, there is a substantial loading of fuel at 
the initial stage unlike the results of the present computations. 

Rext we have obtained complete solutions by varying 
the temperat'ure of the oxidizer stream (T ) and the evaporation 
kinetic parameter (€ ) to estimate their effect on combustion 
process. The general behaviour of the detailed profiles for 
f*, Yq, Yp, ©, and were found to be very similar to those 
obtained for case I, 

The effect of increasing the temperature of the 

external stream (T ) on the interface parameters is shown in 

© 

Pig. 15, Here the data corresponds to case II in Appendix I. 

The surface temperature behaviour is very similar to the one 
observed in case I (see Pig» 12). There occurs a substantial 
increase in the vaporization rate as well as the Interface fuel 
mass fraction compared to the previous case. 

The effect of decreasing the magnitude of evaporation 
constant ( C) on the variables at the interface is shown in 
Pig. 14 corresponding to the data of case III in Appendix I. 










'llllii 
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There occxirs a pronounced reduction in the vaporization rate 
as -well as in the fuel mass fraction at the surface. However, 
the interface temperature remains almost the same in magnitude 
as well as behavio'ur as in case I. At about v = 2 we find 
that vaporization rate has started increasing due to the 
increase in interface temperature in a manner similar to 
case I. However the overall effect is more pronounced here 
because of the fact that there is a reduction in the absolute 
values during the initial time period. Also, it is interesting 
to note that the decreased evaporation constant has also 
reduced the time lag between the building up of with the 

corresponding vaporization rate. 

The combined effect of slow vaporization kinetics 
and increased external stream temperature on the interface 
parameters is shovm in Fig. 15 corresponding to the data of 
case IV in Appendix I, like case III, the slow kinetics has 
reduced vaporization rate substantially as compared to case I 
and thus delaying the building up ef fuel mass fraction inside 


the gaseous boundary layer. Also, as observed earlier higher 
external stream temperature increases the interface temperature. 
This indicates clearly that vaporization kinetics is an 


important factor in determining the vaporization rate. Around 
T = 2 the interface temperature begins to rise more rapidly 
resulting in higher values of the vaporization rate as well as 
fuel mass fraction. It appears that signific|ant-^mio^— 
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activity resiilting in an additional source of heat rises the 
interface temperature after sometime. At around *r = 10 
boiling point of the liquid is reached and the calculations 
continued for the next time inteinrals showed that equilibrixmi 
fuel mass fractions were greater than imity. The present 
analysis then breaks down for further calculations. 

The effect of variation in the relative velocity 
of hot external stream was also studied by varying the 
parameter K. The values of K used in these computations are 
10, 100 4 and 1000, An increase in the value of E actually 
mean's^ an increase in the flow residence time (or low values 
of There was no significant change in the results (see, 

for example, Eig, 13). Computations for case II were also 
done when K = 10 , and here the increase in free stream 
velocity did not show any significant change in the values 
like vaporization rate and interface temperature as compared 
to the resul-ts when K =i 100. But a decreased free stream 
velocity (K = 1000) showed that the fuel loading inside the 
boTondary layer was of smaller magnitude when compared to the 
case when K = 100. This difference in the two cases considered 
can be attributed to the high relative velocity which exists 
between the droplet and the surroxmding hot streaia resulting 
in sweeping out more fuel by convection effects. 

The effect of a choice for initial values of 
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interface temperature (©■_) is significant in the initial 

O 

stages , hut it did not bring any substantial change as far as 
the steady state behaviour is concerned. Similarly the effect 
of initial selection for vaporization rate to start the 
numerical solution was also studied. In Fig, 13 the effect 
of initial guess is also shown. Since no appreciable change 
was foxmd either in the interface temperature or in the 
surface fuel mass fraction development with time or in the 
final steady sta,to values, it was concluded that the initial 
gness can be made somewhat arbitrarily. 

Before concluding the discussion of the results of 
computations it is useful to discuss the time scales which 
could play an important role in the phenomena of combustion 
instability. As mentioned already in Chapter I the unsteady 
interactions can occur if the energy addition is in phase with 
the flow disturbance. In other words, if the product of the 
characteristic time of the combustion process (t_^) and angular 
frequency (co) of the acoustic oscillation is of the order of 
■unity. A typical estimate is made by considering a droplet 
of radius 100 microns injected into a free oxidizing stream 

having a relative velocity of 25 m/sec. Then the characteristic 

-6 • . 

flow time t^ will be of the order of 4 x 10 sec. Making 
use of the time scale at which gaseous phenomena has attained 
steady state (i.e. v =100), the unsteady interactions with 
the gas phase can occur provided that the circular frequency 
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of the aeoustic oscillation is of the same order as tha,t of 

l/t . That is, if « is of the order of 2500 Hz, In general, 
c 

the fuel spray consists of droplets with varying sizes and 
having somewhat different velocities. Consequently, several 
other values of a,coustic frequencies ma,y come into pla.y when 
the comhustion insta^hility problem is analyzed in detail for 
B. rocket motor chanher. 

IV.B CONCLUSIONS 

Detailed profiles for velocity, species mass 
fractions, and temperature ha.ve been obtained for pentaborane 
droplets. In order to investigate the transient b^urning of a 
droplet , it is essential to take into account non-equilibrium 
vaporization kinetics since the steady state situation is not 
obtained for the thermal behaviour within the liquid droplet 
while the gaseous phase phenomena becomes steady much earlier. 
However, the analysis for a droplet burning can be treated as 
a quasi-steady phenomena after a time greater than an order 
of magnitude of the characteristic flow time. 
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APPSBEDIX I 

KljMERIGAI. DATA USED IE COMPUTATIONS 


The physical values for the liq.td.d pentahorane fuel 
used in the numerical computations for different cases are 
summarized below (28-30). 




* 


I 

600 

44.841 

3.0 

540 

o 

H 

100 

and 

1000 

II 

, 1200 

65.245 

3.0 

1080 

10, 

100 

and 

1000 

III 

600 

44.841 

0.1 

540 


100 



IV 

1200 

65.245 

H 

• 

O 

1080 


100 




Cz = 36.12 cal mole"^ deg"^ 

B = 15»000 cal/mole 

H£ = 7200 cal/mole 

q,^ = 1085000 cal/mole 

= 1.987 cal/mole-°K 

T „ = 282°K 

1,0 

The non-dimensional quantities C, Pr , and Sc were 
simply set equa,l to unity. A step size equal to 0.05 was 

found adequate for integration in IQ— direction. To start the 

calcul8.tions in each case we chose the initial va,lues of 
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(Y- ^)° and (-f^) equal to 0.33, 0.1637, 0.33 and 0.1637 
eq ® 

and -3.6, -5.7, -0.116 and -0.183 respectively. 
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